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We report the first observation of two resonance-like structures in the iv + Xci invariant mass 
distribution near 4.1GeV/c 2 in exclusive B° — > K~-K + xa decays. From a Dalitz plot analy- 
sis in which the tv + Xci mass structures are represented by Breit-Wigner resonance amplitudes, 
we determine masses and widths of: Mi = (4051 ± 14±|? ) MeV/c 2 , I\ = (82l^tg) MeV, 
M 2 = (4248l29l 1 !s)MeV/e 2 , and T 2 = (177±||± 3 ^) MeV; and product branching fractions of 
B(B° -> K~Z+ 2 ) x B(Z+ 2 -> tt+Xci) = (3.0±S;l±f:^) x 1(T 5 and (4.0±g;|± 1 g^) x 1(T 5 respectively. 
Here the first uncertainty is statistical, the second is systematic. The significance of each of the 
n + Xci structures exceeds 5 a, including the systematic uncertainty from various fit models. This 
analysis is based on 657 x 10 6 BB events collected at the T(AS) resonance with the Belle detector 
at the KEKB asymmetric-energy e + e~ collider. 



PACS numbers: 14.40.Gx, 12.39.Mk, 13.25.Hw 

INTRODUCTION 

Recently the Belle Collaboration reported the observa- 
tion of a relatively narrow resonance-like structure in the 
ir + ip(2S) mass spectrum produced in B° — > K~ir + yj(2S) 
decays, calling this structure the Z(4430) + [l|. If the 
Z(4430)+ is interpreted as a meson state, then its mini- 
mal quark content must be the exotic combination \ccud) . 
The Z(4430) + observation motivated studies of other 
B° — ► K~ir + {cc) decays. 

In this paper we present a study of the decay B° — * 
K~Tt + Xci, including the first observation of a doubly 
peaked structure in the tt + Xc1 invariant mass distribution 
near 4.1 GeV/c 2 . If the two peaks are meson states, their 
minimal quark content must be the same as that of the 
Z(4430) + . The analysis is performed using data collected 
with the Belle detector at the KEKB asymmetric-energy 
e + e~ collider The data sample consists of 605 fb -1 
accumulated at the T(AS) resonance, which corresponds 
to 657 x 10 6 BB pairs. 



BELLE DETECTOR 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of 
aerogel threshold Cherenkov counters (ACC), a barrel- 
like arrangement of time-of-fiight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) com- 
prising CsI(Tl) crystals located inside a superconducting 
solenoid coil that provides a 1.5 T magnetic field. An iron 



flux-return located outside of the coil is instrumented to 
detect Kl mesons and to identify muons (KLM). The 
detector is described in detail elsewhere |3J. Two dif- 
ferent inner detector configurations were used, a 2.0 cm 
radius beam-pipe and a 3-layer silicon vertex detector for 
the first 155 fb -1 , and a 1.5 cm radius beam- pipe with a 
4-layer vertex detector for the remaining 450 fb -1 

We use a GEANT-bascd Monte Carlo (MC) simula- 
tion [B[ to model the response of the detector, identify po- 
tential backgrounds and determine the acceptance. The 
MC simulation includes run-dependent detector perfor- 
mance variations and background conditions. Signal MC 
events are generated in proportion to the relative lumi- 
nosities of the different running periods. 



EVENT SELECTION 

We select events of the type B° — > K~ir + Xci, where 
the Xci meson is reconstructed via its decay to J/i/j"f, 
with a subsequent J/ip decay to £ + £~ (l + £~ = e + e~ 
or /.t + /i _ ). The inclusion of charge-conjugate modes is 
implied throughout this paper. 

All tracks are required to originate from the beam- 
beam interaction region: dr < 0.2 cm and dz < 2 cm, 
where dr is the distance of closest approach to the beam- 
beam interaction point in the plane perpendicular to the 
beam axis and dz is the corresponding distance along the 
beam direction. Charged pions and kaons are identified 
using a likelihood ratio method that combines informa- 
tion from the TOF system and ACC counters with energy 
loss (dE/dx) measurements from the CDC. The identifi- 
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cation requirements for kaons have an efficiency of 90% 
and a pion misidentification probability of 10%. Muons 
are identified by their range and transverse scattering 
in the KLM. Electrons are identified by the presence of 
a matching ECL cluster with transverse energy profile 
consistent with an electromagnetic shower. In addition, 
charged pions and kaons that are also positively identified 
as electrons are rejected. 

Photons are identified as energy clusters in the ECL 
that have no associated charged tracks detected in the 
CDC, and a shower shape that is consistent with that of 
a photon. 

For J/tp — > e + e~ candidates, photons that have labo- 
ratory frame energies greater than 30 MeV and are within 
50 mrad of the direction of the e + or e~ tracks are 
included in the invariant mass calculation; we require 
\M(e + e-)-m J/i; \ < 50MeV/c 2 . For J/ip -> ^ \x~ can- 
didates we require \M(fi + n~) — mj/^\ < 30MeV/c 2 . To 
enhance the precision of the J/ip energy and momentum 
determination, we perform a mass constrained fit to the 
J/tp candidates. 

For Xci — > J/ipl candidates, we use photons with lab- 
oratory frame energies greater than 50 MeV and require 
\M(J/ipj) - m Xcl | < 30MeV/c 2 . To improve the accu- 
racy of the Xci energy and momentum determination, we 
perform a mass constrained fit to the \ci candidates. 

Candidate B° — * K~n + Xci decays are identified by 
their center-of-mass (cm.) energy difference, AE = 
Y>iEi — -Ebeamj an£ l their beam-energy constrained mass, 
M bc = v^boam - (^Pi) 2 i where E hoam = -Js/2 is the 
beam energy in the cm. and pi and Ei are the three- 
momenta and energies of the B candidate's decay prod- 
ucts. We accept B candidates with 5275MeV/c 2 < 
M bc < 5287MeV/c 2 and \AE\ < 12 MeV. The AE side- 
bands are defined as 24 MeV < |AJ5| < 96 MeV. To have 
well defined Dalitz plot boundaries for both signal and 
sideband events, we perform a mass constrained fit to 
the B° candidates from both regions (to the nominal B° 
mass in all cases). 
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FIG. 1: The AE distribution for the selected B° meson can- 
didates (histogram) and for the \ c \ sidebands (points with 
error bars). The vertical lines indicate the AE signal and 
sideband regions. 
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FIG. 2: The M(J/ipy) distribution for the selected B° meson 
candidates (histogram) and for the AE sidebands (points with 
error bars). The vertical lines indicate the Xci signal and 
sideband regions. 



ANALYSIS OF B° -> K'tt+Xci DECAYS 

The AE distribution for selected B° — > K~ir + Xci can- 
didates is shown in Fig. [TJ The contribution of the Xci 
sideband regions defined as 140MeV/c 2 < \M{J/ip~f) — 
m Xcl \ < 230MeV/c 2 is also shown. The Xci sidebands 
account for almost all the background, which indicates 
that the background is primarily due to combinatorial 
photons; the contamination from events with misiden- 
tified particles is found to be negligibly small. The 
M(J/ip~/) distributions before the Xci mass constrained 
fit for the AE signal and sideband regions are shown in 
Fig. [2] There is a small Xci signal in the AE sidebands 
due to inclusive Xci production in B decays. The J/tp 
signals in the M(/i + /i~) and M(e + e~) distributions are 



almost background-free. 

A signal yield of 2126±56±42 B° ->■ ^"tt+XcI events 
is determined from a fit to the AE distribution using a 
Gaussian function to represent the signal plus a second- 
order polynomial to represent the background. The fitted 
AE resolution, a = (5.93 ± 0.15 ± 0.13) McV/c 2 , is con- 
sistent with the MC expectation of a = (5.62 ± 0.03 ± 
0.09) MeV/c 2 . Here and elsewhere in this report the first 
uncertainty is statistical, the second is systematic. The 
systematic uncertainties for the signal yield and the AE 
width arc estimated by varying the AE interval covered 
by the fit. 

To determine the detection efficiency, we simulate 
B°B° events where B° — > K~ir + Xci with a uniform 
phase-space distribution and the accompanying B° de- 
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cays generically. These MC events are then weighted 
according to the results of the fit to the Dalitz plot that 
is described below. In this way, the reconstruction effi- 
ciency is found to be (20.0 ± 1.4)%, where the following 
sources are included in the uncertainty: the dependence 
on the Dalitz plot model (0.2%); data and MC differences 
for track and 7 reconstruction (1% x 4 for four tracks and 
1.5% for 7), and particle identification (4% for the K~ir + 
pair and 4.2% for £ + £~); uncertainties in the angular 
distributions for Xci — * J/i(H and J/tp — > £ + £~ decays 
(0.2%); and MC statistics (0.6%). The uncertainties from 
different sources are added in quadrature. The efficiency 
is corrected for the difference in lepton identification per- 
formance in data compared to MC, (— 4.5±4.2)%, as esti- 
mated from J/ip — > £ + £~ and e + e _ — > e + e~£ + £~ control 
samples. 

Using (656.7 ± 8.9) x 10 6 as the number of BB pairs 
and Particle Data Group (PDG) 2006 values for the 
branching fractions B(xd —> J/ipj) = 0.356 ± 0.019 and 
B(J/ip^>£+£-) = 0.1187 ±0.0012 0, we determine 

B(B° -> A'-tt+Xci) = (3.83 ± 0.10 ± 0.39) x 10~ 4 . 

The systematic uncertainty includes contributions from 
the uncertainty in the efficiency (7.2%), the systematic 
uncertainty in the signal yield (2.0%), the uncertainty 
due to the variation in the selection requirements (3.9%), 
the uncertainty in the AE signal shape (1.0%, consider- 
ing a sum of two Gaussian functions instead of a sin- 
gle one) and the uncertainties in the Xci and J/ijj decay 
branching fractions (5.3% and 1.0%, respectively). 

The B° -> K-tt + Xc1 decay Dalitz plot (M 2 (ir + X ci) 
versus M 2 (K~ir + )) for the AE signal region is shown 
in Fig. [3] (a). The Dalitz plot distribution exhibits 
some distinct features: a vertical band at M 2 (K~tt + ) ~ 
0.8GeV 2 /c 4 that corresponds to B° -> if*(892)x c i de- 
cays; a clustering of events at M 2 (K~n + ) ~ 2GeV 2 /c 4 
that corresponds primarily to B° — ► if*(1430)x c i decays; 
a distinct horizontal band at M 2 (ir + x c i ) ~ 17GeV 2 /c 4 
corresponding to a structure in the tt + Xc1 channel, de- 
noted by Z + . This latter feature is the subject of this 
report. 

In contrast, the Dalitz plot for the AE sidebands, 
shown in Fig. [3] (b), is relatively smooth and feature- 
less. The Dalitz plot for the phase-space MC candidate 
events, shown in Fig. 01 also exhibits a smooth and fea- 
tureless behaviour. There is a decrease in efficiency in 
the top (bottom) region where the K~ (7r + ) is slow and 
has a low detection efficiency. 

FORMALISM OF DALITZ ANALYSIS 

The decay B° — > K~tt + x c i with the Xci reconstructed 
in the J/ ^7 decay mode and the J/ i/j reconstructed in the 
£ + £~ decay mode is described by six variables (assuming 



23 




M 2 (K V), GeV 2 /c 4 

23 r -o -„ 




M 2 (K V), GeV 2 /c 4 



FIG. 3: The B° -> K'n + x c i decay Dalitz plot for the AE 
signal (a) and sideband (b) regions. 



the widths of the Xci and J/ ip to be negligible) . We take 
these to be M(n+Xci), M(K~-k + ), the xa and J/ip he- 
licity angles {9 Xcl and 0j/^), and the angle between the 
Xci {J/i^) production and decay planes 4> Xcl (</>j/^). Here 
we analyze the B° — > K~n + Xci decay process after in- 
tegrating over the angular variables Xal , Oj/^p, 4>xa an d 
4>j/ip- We find that the reconstruction efficiency is al- 
most uniform over the full 4> Xcl and 4>j/tp angular ranges; 
therefore, after integrating over these angles the inter- 
ference terms between different Xci hclicity states, which 
contain factors of sin </> Xcl , cos </> Xel , sin 2(j> Xal or cos 2(j) Xal , 
are negligibly small. We subsequently verify that the 9 Xcl 
and Oj/if, distributions agree with these expectations. 
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FIG. 4: Dalitz plot for reconstructed B° -> K~tt + Xci MC 
events, generated according to the phase space model. 



We perform a binned likelihood fit to the Dalitz plot 
distribution, where the bin size is chosen by decreasing its 
area until the fit results are unaffected by further changes. 
The selected number of bins is 400 x 400. We consider 
only those bins that are fully contained within the Dalitz 
plot boundaries; this corresponds to 99.3% of the total 
Dalitz plot area. 

In 1.9% of events from the AE signal region we find 
more than one B° candidate. Multiple candidates are 
uniformly distributed over the entire Dalitz plot area. 
No best candidate selection is applied. 

We use a fitting function of the form 

F(s x , s y ) = S(s x , s y ) x e(s x , s y ) + B(s x ,s y ), (1) 

where s x = M 2 (K~-k + ), s y = M 2 {ir + Xci), S and B are 
the signal and background event density functions, and 
e is the detection efficiency. The background B(s x , s y ) is 
determined from the AE sidebands. Its normalization is 
allowed to float in the fit within its corresponding uncer- 
tainty. The bin-by-bin efficiency e(s x ,s y ) is determined 
from the MC simulation. Both sidebands and efficiency 
distributions are smoothed. 

The amplitude for the three-body decay B° — > 
K~ir + Xci is represented as the sum of Brcit-Wigncr 
contributions for different intermediate two-body states. 
This type of description, which is widely used in high 
energy physics for Dalitz plot analyses Q , cannot be ex- 
act since it is neither unitary nor analytic and does not 
take into account a complete description of final state in- 
teractions. Nevertheless, the sum of Breit-Wigner terms 
reflects the main features of the amplitude's behaviour 
and provides a way to find and distinguish the contribu- 
tions of the two-body intermediate states, their mutual 
interference, and their effective resonance parameters. 



Our default fit model includes all known K ir + res- 
onances below 1900MeV/c 2 (re, #*(892), if* (1410), 
#5(1430), K*{U30), #*(1680), # 3 *(1780)) and a sin- 
gle exotic Xci 7r+ resonance. The amplitude for B° — > 
K~tt + Xc1 via a two-body intermediate resonance R (R 
denotes either a K~ir + or n + Xd resonance) and the Xci 
meson in helicity state A is given by 
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Here F { b Lb) and F r Lr 



are the B° meson and R resonance 
decay form factors (the superscript denotes the orbital 
angular momentum of the decay); Mr is the resonance 
mass, sr is the four-momentum-squared and T(s^) is 
the energy-dependent width of the R resonance; T\ is 
the angle-dependent term; (-}^) Lb ■ (~y§^) LR i s a factor 
related to the momentum dependence of the wave func- 
tion, pb (pr) is the B° meson (R resonance) daughter's 
momentum in the B (R) rest frame; and raj is the B° 
meson mass. 

We use the Blatt-Weisskopf form factors given in 
Ref. : 

F (o) =1 
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(2) _ \/ z l + 3z + 9 

Vz 2 + 3z + 9 : 



(3) 



(3) = yjzl + 6z 2 + 45z + 225 
Vz 3 + 6z 2 + 45z + 225 ' 

Here z = r 2 p 2 R where r is the hadron scale, taken to be 
r = 1.6 (GeV/c) -1 , and z = r 2 p 2 R0 where pro is the R 
resonance daughter's momentum calculated for the pole 
mass of the R resonance. For K* resonances with non- 
zero spin, the B decay orbital angular momentum Lb 
can take several values (S, P & D-waves for J = 1; P, D 
& F- waves for J — 2; and D, F & G- waves for J = 3). 
We take the lowest Lg as the default value and consider 
the other possibilities in the systematic uncertainty. The 
energy-dependent width is parameterized as 

T(s R ) = T ■ (pr/ P r ) 2Lr+1 ■ (m R /v/ii) • Fr. (4) 

The angular function T\ is obtained using the helicity 
formalism. For the B° — ► K*(— > K~ ir + )xd decay 

(5) 



Tx = di (e K *), 



where J is the spin of the K* resonance; 6k* 
is the helicity angle of the K* decay. For the 
B° -> K-Z+{-^ tt+Xci) decay 



(6) 



where J is the spin of the Z + resonance and 6 Z + is the 
hclicity angle of the Z + decay. 

In the decays B a -> !£"*(-> A"-7r+)x c i and B° -> 
K~Z + (-^> n + Xci) the parent particles of the Xci are dif- 
ferent and, therefore, the relevant Xci helicity is defined 
relative to different axes: for B° — ► K*{— > A~7r + )x c i 
the axis is parallel to the A~7r + momentum in the Xci 
rest frame; for B° — > K~Z + (-^ tt + Xci) the axis is par- 
allel to the 7r + momentum in the Xci rest frame. The 
angle 9 between the two axes depends upon the event's 
location in the Dalitz plot as indicated in Fig. [5] As a 




M 2 (KV), GeV 2 /c 4 

FIG. 5: cos (9 as a function of the Dalitz plot variables, where 
6 is the angle between the helicity axes for Kir and tt^ci 
intermediate resonances. 

result, the state \ \) z+ with Xci helicity A produced in the 
decay B° — * K~Z + {-^ n + Xci) is not equal to the state 
|A)if» with the same Xci hclicity A produced in the decay 
B° — ► K*(— ► K~n + )xci- The two states are related by 
the Wigner (^-functions via 

|A>*.= X) 4aWIa% + ; (7) 

A'=-1,0,1 

the same relation holds for the amplitudes. 

The resulting expression for the signal event density 
function is 

S(S X , Sy) = 

J2 \j2are^r A f(s x ,s y )+ 

\=-l,0,l K* (8) 

dkW4 + e^' + Af?(s x , Sy )\ 2 , 

A' = -1,0,1 

where and are the normalizations and phases of 
the amplitudes for the intermediate resonance R and Xci 
helicity A. The phase <f>^ ^ 892 ^ is fixed to zero. The 
detector resolution (er ~ 2MeV/c 2 ) is small compared to 
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the width of any of the resonances that are considered 
and is ignored. 

The masses and widths of the K* resonances are fixed 
to their PDG values, except for the k, for which the mass 
and width are allowed to vary within their experimental 
uncertainties Q. The mass and width of the Z + is al- 
lowed to vary without any restrictions. 

RESULTS 

To display the results of the fit, we divide the Dalitz 
plot into four vertical and three horizontal slices as shown 
in Fig.[6j Projections of the fit results for the seven slices 
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FIG. 6: The B° -» K~n + X ci decay Dalitz plot. The solid 
(dashed) lines delimit the four vertical (three horizontal) slices 
that are used to present the fit results in subsequent fig- 
ures. The coordinates of the vertical (horizontal) lines are 
M 2 (K~tt + ) = 1.00GeV 2 /c 4 , 1.75GeV 2 /c 4 and 2.37GeV 2 /c 4 
(M 2 (n + X ci) = 16.0GeV 2 /c 4 and 19.3 GeV 2 /c 4 ). 

are shown in Fig. [3 where the influence of the structure in 
the tt + Xc1 channel is most clearly seen in the second ver- 
tical slice. The mass and width of the Z + found from the 
fit are M = (4150^) McV/c 2 and T = (352+g) MeV; 
the fit fraction of Z + events, defined as the integral of the 
Z + contribution over the Dalitz plot divided by the inte- 
gral of the signal function, S \ A f d ^ d T y \ is (33.lt^)%. 
All quoted uncertainties are statistical. 

The fit fraction is not determined directly from the 
fit and its statistical uncertainty is difficult to estimate 
based on the statistical uncertainties of fit parameters. In 
this paper the statistical uncertainties of fit fractions are 
determined using 1000 toy Monte Carlo samples. Each 
sample is generated according to the probability distri- 
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FIG. 7: Dalitz plot fit with the default model, including one 
Z + . Projections for the slices defined in Fig. [6] are shown: 
(a)-(d) correspond to vertical slices from left to right, (e)- 
(g) correspond to horizontal slices from bottom to top; in (e) 
and (g), plots including the full vertical scale are shown inset. 
The dots with error bars represent data, the solid histograms 
are the fit results, the dashed histograms represent the back- 
ground, and the dotted histograms in (a)-(d) represent the 
sum of all fit components except the Z + . The fitting model 
includes all known K* resonances and one Z+ term. 



bution determined from the fit to experimental data, and 
contains the same number of events as the data. We gen- 
erate 1000 such samples, fit them, and determine the fit 
fractions. We fit the distribution of the obtained fit frac- 
tions to an asymmetric Gaussian function and consider 
the sigmas of the Gaussian function as the statistical un- 



certainty in the fit fraction. 

The significance of the Z + , calculated from the dif- 
ference in log likelihood between fits with and without 
the Z + contribution with the change in the number of 
degrees of freedom taken into account, is 10.7 a. The re- 
sults of the fit with the Z + contribution excluded from 
the default fit function are presented in Fig. [8] The fit 




14 16 18 20 22 



M 2 (x c1 tc + ), GeV 2 



14 16 18 20 22 



M 2 ( Xc X), GeV 2 




M 2 (KV), GeV 2 /c' 



FIG. 8: Dalitz plot fit without any Z + term. Projections 
are shown in (a)-(g) as described for Fig. [7] The dots with 
error bars represent data, the solid histograms are the fit re- 
sults, and the dashed histograms represent the background. 
Here the fitting model with all known K* resonances and no 
resonance in the n + Xci channel is used. 

fractions and significances for each of the resonances in- 
cluded in the default model are listed in Table HI 

To study the model dependence, we consider a variety 



8 



TABLE I: The fit fractions and significances of all contributions for the fit models with a default set of K n + resonances and 
one Z + or two Z + resonances. 





One Z + 


Two Z + 


Contribution 


Fit fraction 


Significance 


Fit fraction 


Significance 


7+ 

z w 


(33.lj*S)% 


10.7a 


(8.0±1;|)% 


5.7(7 


%2 






/in a +6. 1 \ (V 

(10.4I 23 )% 


5.7cr 


K 


(1.9 ±1.8)% 


2.1(7 


(3.6 ± 2.6)% 


3.5(7 


K*(892) 


(28.5 ±2.1)% 


10.6(7 


(30.1 ±2.3)% 


9.8 a 


K* (1410) 


(3.6 ±4.4)% 


1.3(7 


(4.4 ±4.3)% 


2.0 a 


#5(1430) 


(22.4 ± 5.8)% 


3.4cr 


(18.6 ± 5.0)% 


4.5 a 


#1(1430) 


(8.4 ±2.7)% 


5.2 a 


(6.1 ±2.9)% 


5.4 a 


A'* (1680) 


(5.2 ±3.7)% 


2.2 a 


(4.4 ±3.1)% 


2.4 a 


A 3 *(1780) 


(7.4 ± 3.0)% 


3.6 a 


(7.2 ±2.9)% 


3.8(7 



of other fit hypotheses. The results are summarized in 
Table [TTJ The first row in Table [TT| corresponds to the fit 
model with the default set of K~ir + resonances. Rows 2 
through 6 indicate the results from models in which one 
of the Kir resonances: k, K*(U10), #£(1430), If* (1680) 
or #3 (1780), respectively, is removed. Row 7 shows the 
results when a non-resonant Xc\K~ amplitude, param- 
eterized as ae' lb e~ cM ( XclK \ where a, b and c are free 
parameters, is added to the fit model. This amplitude 
can be related to a decay that proceeds via a virtual B* . 
Rows 8 through 10 show results from fits that include the 
non-resonant contribution, but without the #*(1410), 
#*(1680) or #3(1780), respectively. Row 11 corresponds 
to a fit that includes the non-resonant term and releases 
the experimental constraints on the mass and width of 
the k. We also consider models that include the non- 
resonant contribution, plus an additional J = 1 (row 12) 
or J = 2 (row 13) K* resonance with floating mass and 
width. Finally, we replaced the k and A"q(1430) con- 
tributions with the S-w&ve component parameterization 
suggested by the LASS experiment Q (row 14). We used 
the following form of the LASS parameterization (Io| : 

. _ p(l) PB ( >fs 

TTlQ \p (COt — I) 

p 2i5 m l T o/PO \ 

Here s is the four-momentum-squared of the K~tt + pair, 
p is the K~ momentum in the A"~7r + cm. frame, mo 
is the mass and Tq is the width of the AT,* (1430), po is 
the K~ momentum calculated for the pole mass of the 
A"q(1430), and the phase S is determined from the equa- 
tion cot S = ^ + -^ ; where a, b are the model parameters. 
We used the LASS optimal values for the a and b [Toj ] . 

For each fit model the Z + significance is estimated. 
The minimal significance of 6.2 a corresponds to fit model 
13 and is considered as the Z + significance with system- 
atics taken into account. The fit result for model 13 



without the contribution of the Z + is shown in Fig. [S] 
For models with an additional J = 1 or J = 2 AT* 
resonance with floating mass and width, the resulting 
fitted masses and widths of the additional AT* reso- 
nances (M = 2.14GcV/c 2 , T = 3.0 GcV for J = 1 and 
M = 1.05 GeV/c 2 , T = 0.26 GeV for J = 2) do not match 
those of any known Kit resonance Q . 

In the fits described above, the spin of the Z + is as- 
sumed to be 0. We find that the J = 1 assumption does 
not significantly improve the fit quality (in the default fit 
model, —2 In A changes from 17640.7 to 17638.3 for four 
additional degrees of freedom). 

It is not possible to distinguish the contributions of Xci 
hclicity ±1 and —1 in models where the spin of Z + is zero. 
Fits that include both have nearly the same likelihood 
value as fits with only one. 

To address the question of fit quality we constructed 
a two-dimensional histogram with varying bin sizes, in 
which there are a minimum of 16 expected events in each 
bin (95 bins in total). A x 2 is determined for this his- 
togram, x 2 = J2i( n * - fi) 2 / fi-, where is the number 
of events and fi is the expectation value for the i-th bin, 
and a toy MC is used to determine its confidence level. 
For the fit model with the default set of the AT - 7r + res- 
onances and one Z + resonance (Fig. [7]) the confidence 
level is 0.5%. Such a low confidence level value indicates 
that the shape of the structure is not well reproduced 
by a single Brcit-Wigncr. (The confidence levels of the 
fits without a Z + resonance, shown in Figs [5] and [HI are 
3 x 10~ 10 and 9 x 10 -4 , respectively.) 

TWO Z +, S? 

In the Dalitz plot projections for the first and sec- 
ond vertical slices (c/. the top two panels of Fig. [7]) the 
M(xcin + ) structure has a doubly peaked shape. This 
motivated us to add a second Z + resonance to the de- 
fault fit model. The results of the fit with this model are 
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TABLE II: Different fit models that are used to study systematic uncertainties and the significances of the single- and double-^"*" 
hypotheses. 





Model 








Significance 
of one Z + 


One Z + vs. 

r-rA- 

two Z + 


Significance 
of two Z + 


1 


default (see text) 








10.7a 


5.7a 


13.2 a 


2 


no k 








15.6 a 


5.0 a 


16.6 a 


3 


no A" (1410) 








13.4 a 


5.4 a 


14.8 a 


4 


no #0(1430) 








10.4 a 


5.2 a 


14.4 a 


5 


no if* (1680) 








13.3 a 


5.6 a 


14.8 a 


6 


no if* (1780) 








12.9 a 


5.6 a 


14.4 a 


7 


add non-resonant 


XdK 


term 




9.0 a 


5.3 a 


10.3 a 


8 


add non-resonant 


XdK 


" term, no if* (1410) 




11.3a 


5.1a 


13.5 a 


9 


add non-resonant 


XdK 


" term, no if* (1680) 




11.4a 


5.3 a 


13.7a 


10 


add non-resonant 


XdK 


" term, no if* (1780) 




10.8 a 


5.4 a 


13.2 a 


11 


add non-resonant 


XdK 


term, release constraints on k mass & width 


9.5a 


5.3 a 


10.7a 


12 


add non-resonant 


XdK 


term, new if* ( J = 


1) 


7.7 a 


5.4 a 


9.2 a 


13 


add non-resonant 


XdK 


term, new if* ( J = 


2) 


6.2a 


5.6 a 


8.1a 


14 


LASS parameterization 


of S-wave component 




12.4 a 


5.3 a 


13.8 a 



presented in Fig. [TU] 

The confidence level for this fit, calculated using the 
method used for the one Z + hypothesis, is 42%. A com- 
parison of the likelihood values for the one- and two- 
Z + fits favors the two-Z + resonances hypothesis over 
the one-Z + resonance scenario at the 5.7 a level. The 
masses and widths of the two Z + resonances found from 
the two-Z + fit are 

Mi = (4051 ± 14±|?) McV/c 2 , 

r 1 = (82l?^g)MeV, 
M 2 = (4248^™)MeV/c 2 , 

T 2 = (177±H± 3 SS) MeV, 

with fit fractions of /i = (8.0±|;|jl|;|)% and f 2 = 
(^O-'i-f.l- 5 ^)^ - The corresponding product branching 
fractions, calculated as B(B° — ► K~tt + Xci) x fi.2, are 

B(B° -> K~Z+) x B(Z+ -> ir+Xd) = 

(3.0±S;|±?:5) x io- 5 , 

B(B° - K-Z+) x B(Z+ -> n + Xcl ) = 

N.o'n:;;- 1 ;:^ x lo \ 

The product branching fractions are comparable to those 
of the Z(4430) + and other charmonium-like states in a 
leading decay mode [3, [|| . The fit fractions and signifi- 
cances for each of the resonances included in the model 
are listed in Table fl] Wc find that the phase differ- 
ence between the two Z + resonances is close to 7r/2: 

<l>z+ - = L7 -a3- 

To estimate systematic errors, we use the models listed 
in Table [TTl with two Z + resonances instead of one, and 



consider the maximum variations of the Z^ and Z^ 
masses, widths and fit fractions for different fit models as 
a systematic uncertainty. These uncertainties are given 
in the first row of Table IIIIl 

The possibility of multiple minima can be an issue for 
complicated fit models with many contributions. In light 
of this we randomly generated the initial values for the 
fit parameters and repeated each fit several times. The 
deepest minimum is selected. (This approach is used also 
for the single-Z + models.) If any secondary minima are 
within |(5(21nL)| < 2 of the selected solution, they are 
included in the systematic uncertainty determination. 

We also study the systematics due to the uncertainty 
in the form factors for the decays. In addition to the 
default value of the r parameter in the Blatt-Weisskopf 
parameterization r = 1.6 GeV - , wc also consider r = 
1.0 GcV -1 and r = 2.0 GcV -1 . The variation of the Z+ 
parameters is negligible. In addition, we change the as- 
sumption about the value of the B° decay orbital angular 
momentum for those cases where several possibilities ex- 
ist, as discussed above. The resulting uncertainties are 
given in the second row of Table IIIIl 

In the phase-space MC, the angular distributions of 
the Xd —* J/4h an d J/i> — * decays are assumed 

to be uniform. To check the sensitivity of our results 
to this assumption, we weight the MC events according 
to the expectations for the x c iwith zero helicity: 1 + 
2 cos 2 9 Xcl cos 2 Oj/^ — cos 2 9j/^ The variation of the 
Z + parameters is found to be negligible. 

We estimate systematic errors associated with the 
event selection by repeating the analysis while loosening 
the selection requirements on M(J/ipj), M^ c and track 
quality until the background level is increased by a factor 
of two, and while tightening them until the background 
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FIG. 9: An alternative Dalitz plot fit without any Z + term. 
Projections are shown in (a)-(g) as described for Figs[7]and[8] 
point- and line-styles match those from Fig. [8] The fit model 
with all known K* resonances, a \aK~ non-resonant contri- 
bution, and a new KX, but without a Z + term, is used. 



FIG. 10: Dalitz plot fit including two Z + terms. Projections 
are shown in (a)-(g) as described for Fig. [7] point- and line- 
styles also match those of that figure. The fit model with all 
known K* resonances and two Z + terms (Z~{ 



level is decreased by a factor of two compared to the nom- 
inal level. The systematic uncertainties estimated in this 
way are given in the third row of Table IIIII 

In the fits described above, the spins of both Z + res- 
onances are assumed to be zero. We find that a J = 1 
assumption for either or both does not significantly im- 
prove the fit quality, as shown in Table IIVI where we 
show results for all four possible combinations of spin 
J = or J = 1 assignment. The variations in the 
and Z£ parameters for the different spin assignments are 



considered as systematic uncertainties and are listed in 
the fourth row in Table HTT1 

To obtain the total systematic uncertainties, the values 
given in Table ITTT1 arc added in quadrature. 

In the extreme case, i.e. model 2 where the k is elim- 
inated, the two-resonance hypothesis is favored over the 
one-resonance hypothesis with a 5.0 a significance. The 
hypothesis with two Z + resonances is favored over the 
hypothesis with no Z + resonances by at least an 8.1 a 
level for all models. 

We cross-check the estimated significances using toy 
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TABLE III: Systematic uncertainties in the Z + mass, width and fit fraction due to fit model, uncertainty in the form factors, 
uncertainty in the xd, J/tp decay angular distributions and variation of selection requirements. 



Mi, MeV/c 2 



Ti, MeV 



fu% 



M 2 , MeV/c 2 



T 2 , MeV 



Fit model 

Formalism 

Selection 

Spin assignment 



+ is 

-18 

+3 
-0 

+ 
-21 

+ 7 
-30 



+ 15 

- 9 

+8 
-0 

+ 14 

- 

+42 
-20 



+4.6 
-3.0 

+0.4 
- 

+2.4 
-1.4 

+8.0 
-2.6 



+27 
-32 

+ 
-10 

+87 

- 

+ 156 

- 10 



+97 
-34 



+87 

- 

+288 

- 50 



+ 18.5 

- 0.7 

+2.5 

- 

+6.1 

- 

+47.6 




TABLE IV: The — 2 In L values and the change in the number 
of degrees of freedom for the fits with different spin assign- 
ments for the Z± and Z^. 



X 10" 



Thc 



Jl 


J 2 


-21nL 


Ad.o.f. 








17599.2 








1 


17594.3 


4 


1 





17597.5 


4 


1 


1 


17590.1 
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MC. We generated three types of toy MC events accord- 
ing to the fit results of the fit model with the default set 
of K~ir + resonances and with either zero, one or two Z + 
resonances (100 samples of each type). We perform the 
fits to these toy MC samples using the same three fit mod- 
els. The results (mean and r.m.s.) for the significance 
of the single Z + resonance, the level at which the two- 
resonance hypothesis is favored over the one-resonance 
hypothesis and the significance of two resonances com- 
pared to the no-resonance hypothesis for all nine combi- 
nations are given in Table [V] We find that the pattern 



TABLE V: A comparison of the zero, one and two Z + reso- 
nance hypothesis for the toy MC samples generated for 0, 1 
and 2 Z + resonances. The corresponding significances seen 
in the data are given for comparison. 



Hypotheses 


Toy MC samples 




compared 





1 


2 


Data 


1 over 


(1.0 ±0.8) cr 


(9.1 ±1.0) a 


(9.4 ±0.9) (7 


10.7(7 


2 over 1 


(2.0 ±1.2) a 


(1.3 ±0.8) a 


(5.4 ±1.0) a 


5.7a 


2 over 


(1.8 ±0.9) a 


(8.8 ± 1.0) cr 


(10.9 ±1.4) a 


13.2 a 



of the significances observed in data is reproduced well 
by the toy MC with two Z + resonances. 



BRANCHING FRACTION OF THE 

B° -> ,K"*(892) x c i DECAY 

From the K* fit fraction from the two-Z + fit 
given in Table HI we determine the branching fraction 



B(B° - JT(892)°Xci) = (l.T3±g;«±g;i|) 
systematic uncertainty is estimated in the same way 
as described above for the Z^ 2 parameters. The re- 
sult is significantly below the current world average 
(3.2 ±0.6) x 10" 4 @|. However, this is the first measure- 
ment of the branching fraction that takes into account 
interference with other decay channels that produce the 
same final state. The fraction of longitudinal polariza- 
tion is found to be (94.7l!'f lg'g)%, which confirms the 
conclusion that the B — > K* (892)y r1 decay is dominated 



by longitudinal polarization |12t |13| . The significances of 
other intermediate K* resonances are below the 5 a level 
when systematic uncertainties from various fit models are 
taken into account. 



ANGULAR DISTRIBUTIONS OF THE X a AND 
J/ip DECAYS 



Angular distributions for \ci — * and J/ip — > 

£ + £~ decays are not used in the Dalitz analysis and 
therefore provide a useful cross-check. For the \ c \ in 
the helicity zero state the expected angular distribu- 
tion for Xci — * J/4h an d J/ip — > £ + £~ decay is P = 
^(1 + 2 cos 2 9 Xcl cos 2 6 j /^i — cos 2 Oj/^), while for the xa 
in the helicity ±1 state the expected angular distribution 
is Pi = Jtj(1 — cos 2 6 Xcl cos 2 9j/^). Here it is assumed 
that different J/ip helicity states do not interfere. We 
integrate the helicity zero and helicity ±1 components 
of the fit function over the Dalitz plot and find the rel- 
ative contributions wq and w±\. The expected angular 
distribution is then P = wqPo ± w±\P±\. 

The cos# Xcl and cos Q.H4, distributions for the entire 



Dalitz plot, are presented in Fig. QTJ for the leftmost 
vertical slice containing the if* (892) signal in Fig. [L2l 
and for the middle horizontal slice dominated by the Z + 
resonances in Fig. 1131 The agreement with predictions 
is good, ft is evident that the different models give very 
similar predictions and these angular distributions are 
not useful for discriminating between them. 
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FIG. 11: cos 6 Xal (cos 6j/^,) distributions in cos Oj/^ (cos 9 Xcl ) bins for the entire Dalitz plot. The dots with error bars are data, 
the solid (dashed) histograms are the predictions of the model with two J = ( J = 1) Z + resonances, the dotted histograms 
are the predictions of the model with no Z + . The cosBj/^ (cos# Xcl ) bins are ( — 1,-0.6), (—0.6,-0.2), (—0.2,0.2), (0.2,0.6) and 
(0.6,1). 
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FIG. 12: cos# Xcl (cos#j/^,) distributions in cos9j/^ (cos# Xcl ) bins for the vertical Dalitz plot slice that contains the if* (892) 
signal. The dots with error bars are data, the solid (dashed) histograms are the predictions of the model with two J = 
(J = 1) Z + resonances, the dotted histograms are the predictions of the model with no Z + . The cosOj/^, (cos# Xcl ) bins are 
(-1,-0.6), (-0.6,-0.2), (-0.2,0.2), (0.2,0.6) and (0.6,1). 



CONCLUSIONS 



resonance amplitudes, the resonance parameters are 



A broad doubly peaked structure is observed in the 
7r + x c i invariant mass distribution in exclusive B° — > 
K~ir + Xci decays. When fitted with two Breit-Wigncr 



Mi = (4051 ± Utl°) MeV/c 2 , 
ri = (82t^)MeV, 

M 2 = (4248 1 44+180 ^ 

r 2 



29 _ ^)MeV/c^ 



(mt^t^MeV, 



13 




FIG. 13: cos# Xcl (cosOj/^,) distributions in cos9j/^ (cos9 Xcl ) bins for the horizontal slice of the Dalitz plot that contains the 
Z±2 signals. The dots with error bars are data, the solid (dashed) histograms are the predictions of the model with two J = 
(J = 1) Z + resonances, the dotted histograms are the predictions of the model with no Z + . The cos8j/^ (cos# Xcl ) bins are 
(-1,-0.6), (-0.6,-0.2), (-0.2,0.2), (0.2,0.6) and (0.6,1). 



with the product branching fractions of 

B (B° -> K~Zt) x B(Z+ -> ir+xci) = 

(3.o±S:|l?:5) x io- 5 , 



B(B a 



K~Z; 



x B(Z+ 



7T + Xcl) = 



(4.0^-3 + 1^) x 10 - 



The invariant mass distribution M(x c i"" + ) for the Dalitz 
plot slice 1.0GeV 2 /c 4 < M 2 {K-tt+) < 1.75GeV 2 /c 4 , 
where the contribution of the structure in the tt + Xc1 
channel is most clearly seen, is shown in Fig. 1141 

Recently Belle observed the first candidate for a 
charmonium-like state with non-zero electric charge, the 
Z(4430)+ 0]. The two resonance-like structures reported 
here represent additional candidate states of similar char- 
acter. The existence of new resonances decaying into 
XcjTT is expected within the framework of the hadro- 
charmonium model 14j . 

In addition, we measure the branching frac- 
tions B(B° -> K-tt+Xcx) = (3.83 ± 0.10 ± 0.39) x 10~ 4 
and B(B° -> K*(892f Xcl ) = (1.73?%£*o%) x lO" 4 . 
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